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ABSTRACT 

The ability of two-equation turbulence models to accurately predict separated flows is ana- 
lyzed from a combined theoretical and computational standpoint. Turbulent flow past a backward 
facing step is chosen as a test case in an effort to resolve the variety of conflicting results that have 
been published during the past decade concerning the performance of two-equation models. It is 
found that the errors in the reported predictions of the K-e model have two major origins: (1) nu- 
merical problems arising from inadequate resolution, and (2) inaccurate predictions for normal 
Reynolds stress differences arising from the use of an isotropic eddy viscosity. Inadequacies in 
near wall modeling play a substantially smaller role. Detailed calculations are presented which 
strongly indicate that the standard K-e model — when modified with an independently calibrated 
anisotropic eddy viscosity — can yield surprisingly good predictions for the backstep problem. 
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1. INTRODUCTION 


Turbulent flow past a backward-facing step has played a central role in benchmarking the per- 
formance of turbulence models for separated flows. During the past decade — beginning with the 

1980/81 Stanford conference on complex turbulent flows 1 — a variety of two-equation turbulence 
models were tested and compared with the experimental data of Kim, Kline and Johnston and 

Eaton and Johnston 3 for the backstep problem. Initial results 1 indicated that the standard K-e mod- 
el, with wall functions, underpredicted the reattachment point by a substantial amount of the order 
of 20-25%. A significant number of studies have been subsequently published using alternative 
forms of the K-e model wherein a variety of conflicting results have been reported. For example, 

Sindir 4 made modifications to account for streamline curvature based on the algebraic stress model 

of Gibson^; some improvements were obtained although the results were somewhat mixed. Chen^ 
performed calculations with a multiple scale K-e model wherein the turbulent kinetic energy K and 
the turbulent dissipation rate e were decomposed into low and high wavenumber parts along the 

lines suggested by Hanjalic et al? Significantly improved results for the reattachment point were 
obtained (the underprediction was reduced to 5%), however, no detailed comparison of the Rey- 
nolds stresses were made. Speziale^ reported comparable improvements for the backstep problem 
based on an anisotropic K-e model which appeared to suggest that the main source of the errors 

could be due to the use of an isotropic eddy-viscosity in the standard K-e model. Later, So et al 9 
presented results which seemed to indicate that the right near wall model could significantly im- 
prove the predictions of the K-e model for the backstep problem. Similar claims were made by 

Kamiadakis et al 10 who argued that the new near wall treatment in their RNG based K-e model 
gave rise to drastically improved predictions for the reattachment point. However, Avva, Kline 
and Ferziger 1 1 presented results which suggested that the large underprediction of the reattachment 
point attributed to the standard K-e model was due mainly to the fact that the previously reported 
computations were under-resolved. When sufficient resolution was used, they found that the actual 
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error was only of the order of 10% — a deficiency which they eliminated by introducing an empir- 
ical correlation for streamline curvature. Considering the need to accurately predict separated tur- 
bulent flows — which can have a wealth of important scientific and engineering applications — it 
is rather unsettling that such a wide range of conflicting claims still permeate the literature. This 
establishes the motivation of the present paper which is to attempt to clarify this issue. 

In this paper, the following questions will be addressed: 

1) Precisely what does the standard K-e model predict for this backstep problem when the numer- 
ics are done properly to insure full resolution of the flowfield and accurate implementation of 
the boundary conditions, 

2) What is the physical source of actual errors in the standard K-e model, and 

3) Can these errors be eliminated within the framework of two-equation turbulence models with- 
out the ad hoc adjustment of constants or the introduction of other arbitrary empiricisms? 

These questions will be addressed by conducting sufficiently resolved computations of the stan- 
dard and anisotropic K-e model — with several different wall boundary conditions — for the back- 

step configuration of Kim et al ? and Eaton and Johnston 3 . By making detailed comparisons with 
the experiments, it will be demonstrated that the standard K-e model yields results that are within 
12% of the experimental data. With the introduction of an anisotropic eddy viscosity, calibrated 
independently without any ad hoc empiricisms, the K-e model is shown to yield an exceptionally 
good prediction for the reattachment point that is within a few percent of the experimental result. 
The physical implications that these results have on the previous work cited above, as well as on 
future applications to more complex separated flows, will be discussed in detail in the sections to 
follow. 
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2. FORMULATION OF THE PHYSICAL PROBLEM 

The problem to be considered is the fully-developed turbulent flow of an incompressible vis- 
cous fluid past a backward-facing step (a schematic is provided in figure 1). Calculations will be 

conducted for the Kim et al? configuration wherein the expansion ratio (step height: outlet channel 
height) £ is 1:3 and the Reynolds number Re = 132,000 based on the inlet centerline mean velocity 
and outlet channel height. The Reynolds averaged Navier-Stokes and continuity equations are 
solved which take the form: 


9 u dv 
dx + dy " 

9 U -du _9 U _ dp (d 2 U 9 2 u\ dz xx dT xy 

dt +U dx +V ty " ”^ + V l 9 ? + 9/J dx & 

du _ dv —dv _ dp fd 2 v d 2 v A *^xy ^yy 

dt +U dx +V ty ~ <ty + A9? 9y 2 J d y 


( 1 ) 

( 2 ) 

(3) 


where, u and v are the mean velocity components in the x andy directions; p is the modified mean 


pressure;! , x andx are the components of the Reynolds stress tensor x i j = u' i u'y, and vis the 
xx xy yy 

kinematic viscosity. In the standard K-z model with isotropic eddy-viscosity the Reynolds stress 

1 9 

tensor takes the form (see Launder and Spalding ) 


9 — 

x. = iKb.-2C.-Sij 
ij 3 ij »*■ e J 


(4) 


where 


_l(dUi d_Uj' 
ij 2 [dxj dx , . 


(5) 


IS 


the mean rate of strain tensor, K= |x fi is the turbulent kinetic energy, e is the turbulent dissi- 


pation rate, u t = (u, v) is the mean velocity vector, and is a dimensionless constant which is 
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taken to be 0.09. An anisotropic K-e model will be considered — namely, the nonlinear K-e mod- 
el of Speziale 8 — wherein the Reynolds stress tensor takes the form 
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v 3 u ^ e v 
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is the frame-indifferent Oldroyd derivative and Cq is a dimensionless constant which takes on the 
value of 1.68 based on a calibration with turbulent channel flow data. The standard K-e model is 
recovered in the limit as C D —> 0. Both the standard and anisotropic K-e model are solved in con- 
junction with modeled transport equations for the turbulent kinetic energy K and turbulent dissipa- 
tion rate e given by 12 
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where \ T = C ^ — is the eddy viscosity, 
0 


m du du 9u dv 
xx dx x y K dy dx yydy 

is the turbulence production, and C e j, C e 2 , o^- and a e are dimensionless constants which, respec- 
tively, are taken to be 1.44, 1.92, 1.0 and 1.3. The Reynolds averaged equations (l)-(7) are solved 
subject to the following boundary conditions: 

(a) inlet profiles for u, K and e are specified five step heights upstream of the step comer (u is taken 
from the experimental data 3 and the corresponding profiles for K and e are computed from the 
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model formulated for channel flow), 

(b) the law of the wall is used at the upper and the lower walls, and 

(c) standard extrapolated outflow conditions are applied thirty step heights downstream of the step 
comer. 

The law of the wall is applied in two different forms. In the standard two-layer form of the law 
of the wall: 


u + = - lny + + 5, 

Ft 


K _ r -l/2 

2 “ » 
U x 


at the first grid pointy away from the wall if y + =y w T /v > 11.6 given that u + = u/u x (u l is the 

shear velocity and k = 0.41 is the von K£rm£n constant): if y + < 11.6, then u, K and e are inter- 
polated to their wall values based on viscous sublayer constraints. 13 A three-layer law of the wall 
where 


y/ + , for y + £5 

-3.05+5 In y + , for5<y + £30 

5.5 + 2.5 In y + , for;y + >30 

is also implemented along the lines of Avva et al } 1 wherein the fact that the normal derivative of 
K vanishes at the wall is made use of along with more elaborate interpolation formulas for K and 

e (see, also, Amano 13 for a detailed discussion). It must be said at the outset that the law of the 
wall does not formally apply to separated turbulent boundary layers. However, since the separation 
point is fixed at the comer of the backstep — and the flowfield is solved iteratively so that the shear 
velocity u t can be updated until it converges — major errors do not appear to result from its use as 
will be demonstrated in the next section. 

The governing equations (l)-(7) are solved using a finite volume method (see Lilley and 
Rhode 14 and Thangam and Hur 15 ). Here we are interested in the steady state solution which is ob- 
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tained by solving the system of algebraic equations by a line relaxation method 14 with the repeated 
application of the tridiagonal matrix solution algorithm. 16 The issue of resolution is crucial for the 
backstep problem. As discussed by Avva et a/., 11 several researchers have reported results that 

were in severe error due to a lack of adequate resolution. Thangam and Hur 15 conducted a careful 
grid refinement study based on this finite volume method for grids containing 150x75 to 400x200 
mesh points (these meshes have variable grid spacing to allow for maximum resolution near the 
step comer and the walls). The conclusion of their study was that a 200x100 mesh yielded results 
that were within 0.3% of the grid independent solution; the use of significantly coarser meshes 
could lead to appreciable errors whereas, for the most part, the improvement to be gained by the 
use of finer meshes would be minimal. All of the computations conducted in this study were per- 
formed using this 200x100 nonuniform mesh. As indicated earlier, the inlet conditions were spec- 
ified 5 step heights upstream of the step corner and the outlet boundary conditions were specified 
30 step heights downstream of the step comer. It is crucial that a sufficient distance downstream 
of the reattachment point be allowed before the outflow conditions are imposed. Many earlier 
computations of the backstep problem were in significant error due to the imposition of fully-de- 
veloped outflow conditions too close to the reattachment point. 

The steady state solution of Eqs. (l)-(7) is obtained by an iterative solution of the discretized 
equations. The computed solution was assumed to have converged to its steady state when the root 

mean square of the average difference between successive iterations was less than 10~ 4 for the mass 

source. 14,1 ^ Approximately 2000 iterations were needed for the convergence of the standard K-e 
model; this corresponds to approximately 30 minutes of CPU time in a partially vectorized mode 
on the CRAY-2S supercomputer using 64-bit precision. The anisotropic K-e model requires ap- 
proximately 30% more CPU time due to the fact that the additional correction terms to the standard 
K-e model have to be evaluated during each iteration. These correction terms are dispersive rather 
than being dissipative — an additional feature that slows convergence. 
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3. DISCUSSION OF THE RESULTS 


First, results will be presented for the standard K-z model using the two-layer law of the wall 
boundary conditions. For this case — as well as the other results to follow — computed results for 
the mean velocity streamlines, the streamwise mean velocity profiles, the streamwise turbulence 

intensity profiles and the turbulence shear stress profiles will be compared with the Kim et air ex- 
perimental data as updated by Eaton and Johnston. 3 In figure 2(a) the computed streamlines are 
shown which indicate reattachment at X/H ~ 6.0 — a result which is approximately a 15% under- 
prediction of the experimental reattachment point of X/H = 7.1. This is in contrast to earlier report- 
ed results which seemed to indicate that the standard K-z model underpredicts the reattachment 
point by 20-25% — an exaggerated underprediction due to the lack of adequate resolution in those 
computations. (A spurious underprediction of the reattachment point can also result from the ap- 
plication of outflow boundary conditions at a distance X/H < 25 from the step comer as some au- 
thors have done). In figure 2(b), the streamwise mean velocity profiles predicted by the standard 
K-z model are compared with the experimental data. Except in the vicinity of the reattachment 
point, the comparisons are fairly good. More serious discrepancies between the model predictions 
and the experimental data occur in the initial part of the recovery zone for the streamwise turbu- 
lence intensity profiles as shown in figure 3(a). However, the model predictions for the turbulence 
shear stress profiles are reasonably good as can be seen from figure 3(b). 

Next, we will examine the effect of the near wall modeling. In figures 4(a)-(b), the mean flow 
streamlines and streamwise mean velocity profiles obtained from the standard K-z model with the 
three layer law of the wall boundary condition is shown. With this more accurate version of the 
law of the wall, the standard K-z model predicts reattachment at X/H ~ 6.25 — a result which is 

extremely close to that reported by Avva et al . 1 1 By integrating the K-z model directly to the wall 
— with an asymptotically consistent low Reynolds number correction (cf.. So et al . 9 and Speziale 

et a/. 17 ) — it is possible to improve the prediction for the reattachment point by approximately 5%. 
Some earlier published work claiming that the near wall modeling could make a larger difference 
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were in all likelihood due to inadequate resolution of the computations (e.g., So et al . 9 used a 
93x66 mesh). Considering the current state of development of near wall modeling — as well as 
the relatively small benefits that it can yield for the backstep problem — the remaining results to 
be shown will be based on the three-layer law of the wall boundary condition. However, for tur- 
bulent flows where the separation point is not fixed, near wall turbulence modeling can play a much 
more important role. 

Now, it will be demonstrated that the use of the anisotropic eddy-viscosity model (6) can yield 
a more significant improvement in the results. The computed mean velocity streamlines obtained 
from the nonlinear K-e model are shown in figure 6(a). They indicate reattachment at X/H * 6.9: 
a result that is within 3% of the mean experimental value of X/H *7.1. This more accurate predic- 
tion of the reattachment point is reflected in the streamwise mean velocity profiles shown in figure 

6(b) which are in better agreement with the experimental data * than those obtained from the stan- 
dard K-e model based on an isotropic eddy viscosity. Some of the earlier reported results for the 

nonlinear K-e model (see Speziale 8 and Speziale and Ngo 18 ) were lower due to insufficient numer- 
ical resolution of the flowfield. The corresponding streamwise turbulence intensity profiles and 
turbulence shear stress profiles are compared with the experimental data in figure 7(a)-(b). On bal- 
ance, the agreement between the model predictions and experimental data is good (it must be re- 
membered that the uncertainty in the turbulence Reynolds stress measurements is of the order of 
10%). The most notable difference between the predictions of the anisotropic and the standard K- 
e model lies in the streamwise turbulence intensity located approximately one step height above 
the comer of the step. Unlike the standard K-e model, the nonlinear K-e model predicts a dramatic 
trough in this region — a result that is consistent with the more recent independent experi- 
ments. 19,20 In fact, it is the more accurate prediction of the normal Reynolds stress difference 
x - x that accounts for the better predictions of the nonlinear K-e model. This can be illustrated 

yy xx r 

directly in turbulent channel flow. In figure 8, the normal Reynolds stress difference predicted by 
the nonlinear K-e model is compared with the experimental data of Laufer. 21 It is clear that the 
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nonlinear K-e model does a reasonably good job of reproducing the experimental data whereas the 


standard K-e model erroneously predicts that t - x = 0. The accurate prediction of this nor- 

yy xx 

mal Reynolds stress difference is important in the backstep problem since the mean flow stream- 


function V|f (where u = -d\\f/dy , and v = 9\|//0jc) is a solution of the transport equation 


18 


r) r) 
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yy xx 
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( 12 ) 


and, hence, the normal Reynolds stress difference x - X contributes directly to the calculation 

yy xx 

of the mean velocity field. This is the case for any flow with streamline curvature since the normal 
Reynolds stress term on the r.h.s. of (12) can be converted to the alternate form 

d 2 (x -x ) 

v nn sa ' 

dnds 

in terms of an intrinsic coordinate system where n and s are, respectively, the normal and tangential 
directions along a streamline. 

Finally, we will examine the RNG based K- 8 model of Yakhot and Orszag. 22 For high turbu- 
lence Reynolds numbers, it is of the same general form as the standard K-e model given by Eqs. 
(4), (8) and (9). However, the numerical values of the constants — which are computed by the 
RNG approach — take on the alternate values: 

C|j= 0.0837, C el = 1.063, C e2 = 1.7215 
= 0.7179, c e = 0.7 179, K = 0.372 

For the low turbulence Reynolds numbers which occur close to the walls in this backstep problem, 
the RNG K-e model has a built-in correction that allows for a direct integration to a solid boundary 
with the no-slip condition applied. However, we found this near wall correction to be rather am- 
biguous — especially in the presence of geometrical discontinuities such as those which occur near 
the backstep comers. Hence, we will present results in which we match to the three-layer law of 


- 9 - 


the wall. In figures 9(a)-(b), the mean flow streamlines and the streamwise mean velocity profiles 
obtained from the RNG K-z model are shown. The RNG K-z model predicts reattachment at 
X/H ~ 4 a result that constitutes a substantial underprediction of the experimental reattachment 
point of X/H ~ 7.1. Due to this significant underprediction of the reattachment point, the stream- 
wise turbulence intensity and turbulence shear stress profiles predicted by the RNG K-z model are 
in serious error as illustrated in figures 10(a)-(b). The origin of this substantial underprediction of 
the reattachment point lies in the choice of C E j = 1.063. In a homogeneous shear flow, it can be 
shown that the eddy viscosity 


v T ~exp (Xt*) 


(13) 


where t* is the time (nondimensionalized by the shear rate) and X is the growth rate given by 23 


X = 


(C £ i ~ 1) (C e2 — 1) 


, 1/2 


(14) 


Hence the growth rate of the eddy viscosity becomes singular as C el 1. For C el = 1.063 as 

opposed to the more standard value of C e j = 1.44 — the growth rate of the eddy viscosity will be 
overly large. An overprediction of the eddy viscosity will make the model too dissipative a fea- 

ture which will cause the separation bubble to shorten as shown in figure 9(a). Hence, while the 
use of an alternate near wall treatment with the RNG K-z model could cause the separation bubble 
to become slightly larger, we feel that it is very unlikely that a reattachment point close to X/H ~ 7 

could be obtained from this RNG K-z model as reported in Kamiadakis et al}® (We believe that 
the calculations of Kamiadakis et cl}® may have been under-resolved). It should be noted, how- 
ever, that Smith and Yakhot 24 have recently determined that there may have been an error in the 
original calculation of C el by Yakhot and Orszag 22 ; they obtained an alternate value of C el = 1.40. 
With this higher value of C e j and with the use of an RNG based anisotropic eddy-viscosity (see 

Rubinstein and Barton 25 ) — there is no question that the RNG K-z model would yield excellent 
results for the backstep problem. 
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4. CONCLUSIONS 


A reasonably detailed theoretical and computational analysis of the performance of two-equa- 
tion turbulence models for turbulent flow past a backward facing step has been presented in an ef- 
fort to resolve numerous conflicting studies that have been published during the past ten years. The 
following general conclusions were arrived at: 

1) The standard K-e model with three-layer wall functions predicts reattachment at X/H = 6.25 — 
a result that is within 12% of the experimental value of Kim, Kline and Johnston 2 for this back- 

step problem. This result is in close agreement with the recent computations of Avva et al . 
and hence there is little doubt that earlier studies which attributed a 20-25% underprediction of 
the reattachment point to the standard K-e model were in error due to inadequate numerical res- 
olution of the flowfield. 

2) When the standard K-e model is appropriately modified to include an anisotropic eddy viscos- 
ity, 8 the reattachment point moves out to X/H « 6.9. This result, as well as the other mean ve- 
locity statistics, are within 3% of the experimental data — an improvement due to the more ac- 
curate prediction of normal Reynolds stress differences. 

3) The RNG K-e model of Yakhot and Orszag 22 substantially underpredicts the reattachment 
point. This problem arises since C e \ is too close to 1 — a feature that causes the model to be 
overly dissipative (e.g., when C el = 1, the growth rate of the eddy viscosity becomes singular 

in homogeneous shear flow). The previously published results of Kamiadakis et al.}® to the 
contrary, appear to have been under-resolved. However, there is no question that with the re- 
cent suggested correction (Smith and Yakhot 2 ^) of Cgj to 1.40 — and with the use of an aniso- 
tropic eddy viscosity — the RNG K-e model would yield excellent results for the backstep 
problem. 

Finally, some remarks are in order concerning the implications that these results have on tur- 
bulence modeling. The deficiencies of two-equation models are well established, particularly in 

turbulent flows with body forces or Reynolds stress relaxation effects. 2 ^ Consequently, the find- 
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ings of this study should not be interpreted as an unequivocal endorsement of two-equation models. 
Nonetheless, this study certainly does indicate that properly calibrated two-equation turbulence 
models, with an anisotropic eddy viscosity, can yield excellent results for the backstep problem that 
are far superior to the results obtained from the older zero or one equation models. 
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List Of Figures 

Figure 1 Physical Configuration and Coordinate System. 

Figure 2 Computed Flowfield for the Standard K-z Model [E = 1:3; Re = 132,000; 200x100 
mesh; C ^ = 0.09; C e j = 1.44; C e 2 = 1.92; a K = 1.0; o e = 1.3; k = 0.41] 

a) Contours of mean streamlines 

b) Mean velocity profiles at selected locations ( computed solutions based on two 

layer wall model; 0 experiments of Kim et al, 1980; Eaton & Johnston, 1981) 

Figure 3 Computed Turbulence Stresses for the Standard K-z Model [E = 1:3; Re = 132,000; 

200x100 mesh; = 0.09; C el = 1.44; C e2 = 1.92; a K = 1.0; c e = 1.3; k = 0.41; 

computed solutions based on two layer wall model; o experiments of Kim et al, 1980; 
Eaton & Johnston, 1981] 

a) Turbulence intensity profiles 

b) Turbulence shear stress profiles 

Figure 4 Computed Flowfield for the Standard K-z Model [E = 1:3; Re = 132,000; 200x100 
mesh; = 0.09; C e i = 1.44; C e2 = 1.92; o K = 1.0; o e = 1.3; k = 0.41] 

a) Contours of mean streamlines 

b) Mean velocity profiles at selected locations ( computed solutions based on three 

layer wall model; 0 experiments of Kim et al, 1980; Eaton & Johnston, 1981) 

Figure 5 Computed Turbulence Stresses for the Standard K-z Model [E = 1:3; Re = 132,000; 

200x100 mesh; = 0.09; C E i = 1.44; C e2 = 1.92; 0 % = 1.0; a e = 1.3; k = 0.41; 

computed solutions based on three layer wall model; o experiments of Kim et al, 1980; 
Eaton & Johnston, 1981] 

a) Turbulence intensity profiles 

b) Turbulence shear stress profiles 

Figure 6 Computed Flowfield for the Nonlinear K-z Model [E = 1:3; Re = 132,000; 200x100 
mesh; = 0.09; C e i = 1.44; C e2 = 1.92; a K = 1.0; a e = 1.3; k = 0.41; C D = C E = 1.68] 

a) Contours of mean streamlines 

b) Mean velocity profiles at selected locations ( computed solutions based on three 

layer wall model; 0 experiments of Kim et al, 1980; Eaton & Johnston, 1981) 
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Figure 7 Computed Turbulence Stresses for the Nonlinear K-e Model [E = 1:3; Re = 132,000; 

200x100 mesh; = 0.09; C e j = 1.44; C e 2 = 1.92; - 1.0; o e = 1.3; k = 0.41; Cq = 

C E = 1.68; computed solutions based on three layer wall model; 0 experiments of 

Kim et al, 1980; Eaton & Johnston, 1981] 

a) Turbulence intensity profiles 

b) Turbulence shear stress profiles 

Figure 8 Comparison of the Computed Values of the Normal Reynolds Stress Difference Ob- 
tained from the Nonlinear K-e Model with the Experimental Data of Laufer (1951) for 
Fully-Developed Turbulent Channel Row. (It should be noted that the standard K-e 
model predicts a zero normal stress difference) 

Figure 9 Computed Flowfield for the RNG K-e Model \E = 1:3; Re = 132,000; 200x100 mesh; 
= 0.0837; C el = 1.063; C e2 = 1.7215; c K = a E = 0.7179; k = 0.372] 

a) Contours of mean streamlines 

b) Mean velocity profiles at selected locations ( computed solutions based on three 

layer wall model; 0 experiments of Kim et al, 1980; Eaton & Johnston, 1981) 

Figure 10 Computed Turbulence Stresses for the RNG K-e Model [E = 1:3; Re - 132,000; 

200x100 mesh; C ^ = 0.0837; C el = 1.063; C e2 = 1.7215; <5 K = a e = 0.7179; K = 0.372; 

computed solutions based on three layer wall model; 0 experiments of Kim et al, 

1980; Eaton & Johnston, 1981] 

a) Turbulence intensity profiles 

b) Turbulence shear stress profiles 
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Physical Configuration and Coordinate System 

Figure 1 
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(a) Streamlines 



U o 


(b) Dimensionless Mean Velocity Profile 

( Computations with the two layer model 

o Experiments of Kim etal, 1980; Eaton & Johnston, 1981) 

Computed Flowfield for the Standard K-e Model 
[E = 1:3; Re = 132,000; 200x100 mesh; C H = 0.09; C el = 1.44; C e2 = 1.92; 
o K = 1.0; a e = 1.3; K = 0.41] 

Figure 2 


- 17 - 


5.33 


7.67 


8.55 


10.33 


15.66 



(a) Turbulence Intensity 



(b) Turbulence Shear Stress 


Computed Turbulence Stresses for the Standard K-e Model 
[E = 1:3; Re = 132,000; 200x100 mesh; = 0.09; C el = 1.44; C e2 = 1.92; 

&K- "I -0> <^e = 1 *3; k = 0.41 ; Computations with the two layer model; 

o Experiments of Kim etal, 1980; Eaton & Johnston, 1981] 

Figure 3 
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Uo 


(b) Dimensionless Mean Velocity Profile 

( Computations with the three layer model; 

o Experiments of Kim et al, 1980; Eaton & Johnston, 1981) 

Computed Flowfield for the Standard K-e Model 
[E = 1:3; Re = 132,000; 200x100 mesh; C^ = 0.09; C el = 1.44; C e2 = 1.92; 
g k = 1.0; o e = 1.3; k = 0.41] 

Figure 4 
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(b) Turbulence Shear Stress 


Computed Turbulence Stresses for the Standard K-e Model [E = 1:3 
Re = 132,000; 200x100 mesh; = 0.09; C el = 1.44; C e2 = 1.92; <j K = 1.0 

G e = 1.3; k = 0.41; Computations with the three layer model 

o Experiments of Kim etal, 1980; Eaton & Johnston, 1981] 

Figure 5 
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(a) Streamlines 


X/H : 1.33 2.67 5.33 6.22 8.00 10.67 



Uo 


(b) Dimensionless Mean Velocity Profile 

( Computations with the three layer model; 

o Experiments of Kim et al, 1980; Eaton & Johnston, 1981) 


Computed Flowfield for the Nonlinear K-e Model 
[E = 1:3; Re = 132,000; 200x100 mesh; C^ = 0.09; C eI = 1.44; C e2 = 1.92; 
(j/e ~ "1 -0; (J £ = 1.3; k = 0.41 ; Cq = Cg = 1 .68] 

Figure 6 
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(a) Turbulence Intensity 


X/H : 1.33 5.33 7.67 8.55 10.33 15.66 



(b) Turbulence Shear Stress 


Computed Turbulence Stresses for the Nonlinear K-e Model [E = 1:3; 
Re = 132,000; 200x100 mesh; = 0.09; C el = 1.44; C e 2 = 1.92; g k = 1.0; 

a e = 1 .3; k = 0.41 ; C D = C E = 1 .68; Computations with the three layer 

model; o Experiments of Kim et al, 1980; Eaton & Johnston, 1981] 

Figure 7 
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Comparison of the computed values of the normal Reynolds stress dif- 
ference obtained from the nonlinear fC-s model with the experimental 
data of Laufer (1951) for the fully-developed turbulent channel flow. (It 
should be noted that the standard K-e model predicts a zero normal 
stress difference) 

Figure 8 
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(b) Dimensionless Mean Velocity Profile 

( Computations with the three layer model; 

o Experiments of Kim etal, 1980; Eaton & Johnston, 1981) 


Computed Flowfield for the RNG K-e Model 
[E = 1:3; Re = 132,000; 200x100 mesh; 

C e2 = 1.7215; O k = 0.7179; G e = 0.7179; K = 0 

Figure 9 


C„ = 0.0837; C el = 1.063; 
372] 
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(a) Turbulence Intensity 


X/H : 1.33 5.33 7.67 8.55 10.33 



(b) Turbulence Shear Stress 

Computed Turbulence Stresses for the RNG K-e Model [E = 1:3; 
Re = 132,000; 200x100 mesh; = 0.0837; C el = 1.063; C e2 = 1.7215; 

o K = 0.7179; a e = 0.7179; k = 0.372; Computations with the three 

layer model; o Experiments of Kim etal, 1980; Eaton & Johnston, 1981] 

Figure 10 
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